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Motivation 
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Combustion  instability  is  an  unsustainable  cimwth 
of  pressure  and_  heat  transfer  fluctuations  in  a 

rocket  engine 

—  Irreparable  damage  can  occur  in  <1s 

Combustion  Instability  caused  a  4-yr  delay  in  the 
development  of  the  F-1  engine  used  in  the  Apollo 
program 

—  More  than  $400M  for  the  propellants  alone  at 
2010  prices 

—  More  than  2000  full  scale  test 


Damaged  engine  injector 
faceplate  caused  by  combustion 
instabilities 


“Combustion  instabilities  have  been  observed  in  almost  every 
engine  development  effort,  including  even  the  most  recent 
development  programs”  -  current  JANNAF  Stability  Panel  Draft 
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Motivation  (cont’d) 
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Statement  of  need 


Modern  rockets  operate  at  supercritical  pressures  with 
respect  to  the  propellants 

*  Need  to  understand  mixing  and  combustion  beyond 
liquid,  gas  states 

Shear  coaxial  injectors  are  a  common  choice  for  cryogenic 
liquid  rocket  engines 

Interactions  of  transverse  acoustics  with  injector’s  own 
modes  and  mixing  needs  to  be  understood  for  combustion 
instability 

*  Need  to  understand  differences  in  response  to  pressure 
and  velocity  nodes 

Understand  what  non-dimensional  numbers  capture  the 
mixing  of  typical  injectors 

Characterize  how  geometry  affects  mixing 
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Highlights  from  previous  work  on  jet 
instabilities  1/2 


Michalke,  1964 

—  Linear  stability  theory  for  inviscid  instability  of  a 
hyperbolic  tangent  velocity  profile 

Crow  and  Champagne,  1971 

—  Single  jet  preferred  mode,  Std=fd/U~0.3 
Ko  etal, 1976-1989 
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—  Some  of  earliest  detailed  description  of  near  field  mixing 
for  coaxial  jets 

Boldman  et  al,  1975 

—  Experimental  and  theoretical  analysis  for  mixing  of  two 
air  streams  with  different  velocities  -  points  out  different 
vortex  interactions,  St,~0.2  (Uave) 

Gutmarkand  Ho,  1983 

—  Collects  previous  results  on  jet  preferred  mode,  Std  has  a 
range  from  ~0.24-0.64 


Kwan  and  Ko,  J.  Sound  and  Vibration,  48  (2),  1976 


Wicker  and  Eaton,  AIAA  J,  (32)  No.3,1994 


•  Wicker  and  Eaton,  1994 


—  Forces  air  inner  and  outer  jets  independently  -  observes 
vortex  growth 

*  Dahm  et  al,  1992 


—  Seminal  pictures  of  different  instabilities  plus  effect  of 
absolute  velocity  and  R 
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Highlights  from  previous  work  2/2 
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Balarac,  da  Silva,  Metais  et  al  (2003,  2007) 

—  DNS  analysis  of  coaxial  jets  -  same  density,  top-hat  profiles 

—  Consider  two  shear  layers,  study  effect  of  R 

^  ^  7777777771 

—  Consider  axisymmetric  and  azimuthal  excitation 
Buresti,  Talamelli,  Petagna  (1994, 1998) 

—  Air  jets,  same  density,  top-hat  profile,  Stdo=fdo/Uoj~0.3  to  1 , 
function  of  xlDt 

Segalini,  Talamelli,  et  al  (2006,  2011) 

—  Air  jets,  same  density,  top-hat  profile,  Stb(Mp)=f b( I i p)/Uaverage 
Birbaud,  Ducruix,  Durox,  Candel  (2006-2007) 

—  Single  air  jets,  low  Re,  laminar,  top-hat  profile,  subjected  to 
acoustic  modulation 

—  Systematic  study  of  effect  of  modulation  in  terms  of  Std,  St0 
Tshohas,  Canino,  Heister  (2004,  2009) 

—  2D  unsteady  CFD  for  LOX/H2  elements  but  non-reacting 

—  Unforced  behavior,  found  StMp=fd|OX/U|OX~0.1 0-0.25 
Richecoeur,  Scouflaire,  Ducruix,  Candel  (2006) 

—  Forced  transverse  acoustic  excitation  of  flames 
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Balarac  et  al,  Phys  of  Fluids  (19),  2007 


Segalini  et 
al,  Phys  of 
Fluids  (23), 
2011 
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Richecoeur  et  al,  JPP  (22)  No  4,  2006 
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Thermodynamic  States  of  N2 
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Property  gradients  around  critical 
point 
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Geometric  Configurations 


LAR-thin 


Similar 

lip 

thickness 


Similar 


SAR-thin 
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SAR-thick 


FsAR,  LAR  Small,  Large  Area 
Ratio 

Thick,  Thin  ->•  Post  lip  thickness 

Davis,  Rodriguez,  Leyva  et  al.  (0.5  D1  recess) 
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Rodriguez  (0.5  D1  recess),  Graham  etal.  (no 
recess) 
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Relevant  variables  for  cold-flow  studies 


Geometry 

Acoustics 

Recess 

Phase 

Single  Jet 

On/Off 

N/A 

v’  max 

LAR_thickLip 
(injector  1) 

1/2D1 

SAR_thinLip 
(injector  2) 

Off 

P’  max 

U’max 

1/2D1 

D1 

0 

2-phase  P<Pc 

SAR_thickLip 
(injector  3) 

-0.2 

P>Pc 

T>Tc 

T<Tc 

LAR_thinLip 
(injector  4) 

-0.2 
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N/A 


0.1-20 

0.1-20 


Coupling 


No 


No 


Yes&  No 
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Experimental  setup  -  EC-4 
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Piezo-Siren 


Coaxial  Injector 


Chamber  Pressure  vs.  Time 
(Subcritical,  J  =  9.6,  f  =3.01  kHz) 


0  0.5  1  1.5  2  2.5  3  3.5 

time  (ms) 


Thermocouple  and  P 
transducer 


Supercritical 
fluid  region 


Dark-Core  Length  Measurement 
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•  First  raw  grayscale  images  were  converted  to  binary  images 

•  A  contour  was  drawn  around  the  “dark-column”  in  the  binary  image 

•  Axial  length  of  the  dark-column  measured  and  defined  as  the  Dark-Core  Length,  L 

Grayscale  (8-bit)  Binary  Image  Raw  Image  w / 
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State  of  the  art  for  mixing  in  shear 
coaxial  jets  -  2007 


O  Davis  Subcritical  P;  High  T  (*) 
X  Davis  Subcritical  P;  Low  T  (*) 
XFavre-Marinet  DR=0.138  Air 

■  Favre-Marinet  DR=0.028  He 
+  Au  and  Ko 

♦  Subcritical  P,  High  T 

■  Nearcritical  P,  Low  T 


O  Davis  Nearcritical  P;  High  T  (*) 
X  Davis  Nearcritical  P;  Low  T  (*) 

■  Favre-Marinet  DR=0.655  Air 
□  Rehab  et  al.  D3/D1  =  1.37 

x  Woodward  KI(aq)  -  N2 

■  Subcritical  P,  Low  T 

♦  Supercritical  P,  High  T 


O  Davis  Supercritical  P;  High  T  (*) 
XEroglu  et  al.  Re=9328 

■  Favre-Marinet  DR=0.138  He 

■  Rehab  et  al.  D3/D1  =2.29 

•  Woodward  KI(aq)  -  He 

•  Nearcritical  P,  High  T 

•  Nearcritical  P,  High  T  (2) 


0.01  0.1  1  10  100  1000 


Outer  to  Inner  Jet  Momentum  Flux  Ratio,  J 
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Dark  Core  Length  -  SAR-thick  -  All 
data 
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Dark  Core  Length  without  Acoustics 


Outer  to  Inner  Jet  Momentum  Ratio  (J) 


■  Inj  1  Sub  Davis 

■  Inj  1  Super,  Davis 

♦  Inj  1  Sub,  Leyva 

♦  Inj  1  Near,  Leyva 

♦  Inj  1  Super,  Leyva 

▲  Inj  1  Sub,  Rodriguez 

▲  Inj  1  Near,  Rodriguez 

▲  Inj  1  Super,  Rodriguez 
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Dark  Core  Length  -  SAR-thick  -  All 
data 
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Dark  Core  Length  without  Acoustics 


■  Inj  1  Sub  Davis 

□  Inj  1  Tr<1  Near,  Davis 

■  Inj  1  Super,  Davis 

♦  Inj  1  Sub,  Leyva 

♦  Inj  1  Near,  Leyva 

o  Inj  1Tr<1  Near,  Leyva 

♦  Inj  1  Super,  Leyva 

A  Inj  1  Sub,  Rodriguez 

a  Inj  1  Near,  Rodriguez 

A  Inj  1  Tr<1  Near, 
Rodriguez 

A  Inj  1  Super,  Rodriguez 


Outer  to  Inner  Jet  Momentum  Ratio  (J) 
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Dark  Core  Length  -  SAR  thick 
supercritical  pressure 


Dark  Core  Length  -  SAR  thick  - 
supercritical  pressure 
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Geometric  Configurations 
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Baseline  Flows:  LAR-thin  Injector 
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Baseline  Flows:  SAR-thick  Injector 
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Sub  and  Nearcritical  Baseline  Dark 
Core  Lengths  for  4  Geometries 
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Baseline  Normalized  DCL,  LBID1  - 
Pr  =  0.44 
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Baseline  Normalized  DCL,  L,  ID  - 
Pr=  1.05 
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Baseline  Normalized  DCL,  LBID1 
Collapse 


•  Non-linear  regression  of  LB!DX  data  from  all  four  injectors  at  each  Pr 


C4 


Pr  =  0.44 


G 


1.05  11  -0.43  -0.12  0.15 
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ACOUSTIC  EXCITATION 

ANALYSIS 
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Acoustic  Field  Set-Up:  Pressure 
Antinode 


*  Pressure  antinode  (PAN)  -  condition  of  maximum  pressure  perturbation 
in  the  acoustic  field 

*  Piezo-sirens  forced  in-phase 

*  Superposition  of  quasi-ID  acoustic  waves  traveling  in  opposite 
directions  =>  PAN  at  the  jet  location  (geometric  center  of  test  section) 


t  (ms) 


t  (ms) 
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Image  interpretation  key 
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Baseline: 

Acoustics 
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-Critical  Pressure:  Two  Geometries 
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ear-Critical  Pressure:  Two  geometries 


LAR thick 


LAR  thin 


SAR_thinLip 
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Pr=1.05,  J=l.l,  p7p=0.32% 
Coupled  to  Outer  Jet  Acoustic  Mode 


Nearcritical  conditions 
*  No  jet  bending 
observed 

•  large  vortical 

structures 
generated  when 
coupled  to  injector 
mode 

*  Reduction  of  dark 
core  can  be  as  large 
as  90% 

•  More  clear  response 
of  jet  to  pressure 
antinode 

Subcritical  conditions 

•  Same  mode  - 
vortical  structures 

*  Not  as  dramatic 
reduction  as  with 
nearcritical  cases 


Pr=1.05,  J=2.2, 
p7p=0.25% 

Non  Coupled 
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Amplitude  Amplitude  Amplitude 


New  analysis:  synchronized  p’  and 
images  taken  with  microscopic  lens 
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Inner  Jet  Pressure 


Chamber  Wall  Center  Pressure 


Outer  Jet  Pressure 


Pr  =  1.03,  J  =  5.1,  VR  =  5.9 
Forcing  Condition:  3000  Hz,  Maximum  Ap 


Uo=6.  25  m/s,  Ui=1.06  m/s,  To=165K  ,  Ti=118  K 
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Results  -  Baseline  Low  J  Flow  at  Pr 
=  0.44 


•  LAR-thin,  Pr=  0.44,  J  = 
0.5 

Snapshot  Average 


POM  1  POM  2 


Antisymmetric 
Structures  Identified  with 
Characteristic  Frequencies 


Amplitude  information 
contained  in  singular  values 


10  10 
Mode 


Power  Spectral 
Densities  (PSD) 
of  Temporal 
Coefficients  of 
POMs  1  and  2 
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Sample  Animation  -  PAN  (fF  =  3.14 
kHz) 


•  LARthin  Pr=  0.44, 
J  =  0.5 


Raw  Images  ROMs  1  and  2  POM  1  POM  2 

Time  =  0.00  ms  Time  =  0.00  ms  Time  =  0.00  ms  Time  =  0.00  ms 
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CPSD  Magnitude 


Results  -  LAR-thin,  Pr  =  0.44,  Baseline  f  © 
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*  Antisymmetric  flow  structures  indicated  helical  type  flow  instabilities  for 

J=  2.1  J=  5.2  ,/=  11  J=  20 


Characteristic  peaks  broadened  and  shifted  to  higher 
frequencies  with  increasing  outer  jet  velocity 


■  oivtic*  I 
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CPSD  Magnitude 


Results  -  LAR-thin,  Pr  =  0.44,  PAN 


Gradual  shift  from  symmetric  to  antisymmetric  flow  structures  with 
increasing  J 


•  Response  at  forcing  frequency,  fF,  dominant  at  lower  J 

j=  2.1  J  =  5.2  /=n  J=  20 


(fv  =  3.12  kHz)  %  =  3.12  kHz)  (fF=  3.10  kHz)  (/F  =  3.11  kHz) 
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Results  -  SAR-thin,  Pr=  0.44 
Baseline 


J 
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9  Helical  type  flow  instabilities  became  more  well-defined  with  increasing  J 


12 


4  D 


8  D 


J=  17 


J  =  2.0 


io6 


/(kHz) 


Unlike  LAR  flows,  characteristic  peaks  showed  minimal  variation  in 
_ frequency  with  outer  jet  velocity _ 
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CPSD  Magnitude 


Results  -  SAR-thin,  Pr=  0.44,  PAN 


*  Symmetric  structures  persist  despite  increasing  J 

•  Response  seen  for  all  J 


J-  2.0  J  =  5.2  J  =  12  J  =  17 

(/f  =  2.97  kHz)  (/f  =  3.02  kHz)  (/F  =  2.92  kHz)  (fF  =  2.90  kHz) 
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Results  -  SAR-thick  -  Pr  =  0.44  Baseline 


*  Well-defined  antisymmetric  flow  structures  observed  for  J  >  2 

•  Unlike  the  LAR-thin  flows,  peaks  in  CPSD  magnitude  spectra  remained  at 
the  lower  end  of  the  frequency  spectrum  despite  increasing  J 
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Results  -  SAR-thick  -  Pr  =  0.44,  PAN 


39 


J=  0.5  J=  2.1  J-  5.7 


J=  10  J  =  21 


Collaboration  with  ECP:  Thomas  Schmitt,  Juan  Rodriguez 
(AFRL  post-doc),  Sebastien  Candel,  Ivett  Leyva 
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Velocity  Fields  for  Different  J  values 


Top:  centerline  profile  of  density;  bottom: 
centerline  profile  of  velocity.  Dark  line:  J  = 
1.1;  dash  point  line:  J  =  3.0;  dashed  line:  J  = 
9.3. 


Longitudinal  slice  of  axial  velocity.  White:  minimum;  black:  maximum. 
Dark  line  indicates  iso-contour  of  zero  axial  velocity.  Top:  J  =  1.1;  middle:  J 

=  3.0;  bottom:  J  =  9.3. 
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Fundamental  frequencies  for 
baseline  conditions 


Case  “N6” 
(J=3.05) 
Density 
distribution 
(white:  60 
kg/m3; 
black:  410 
kg  /m3; 
logarithmic 
scale). 
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0.14 

Case  “N8” 

(J  =  9.3) 
Density 
distribution 
(white:  60  kg/ 
m3;  black: 
maximum ; 
logarithmic 
scale). 
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Found  relevant  St  numbers  for  our  configurations 
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*  The  single  largest  unknown  in  combustion  instability 
is  the  “combustion  response”  (how  combustion 
responds  to  acoustic  waves) 

*  Within  the  combustion  response,  flame  holding  in  the 
near  injector  field  is  a  key  mechanism 

*  GA  Tech  was  selected  to  lead  an  effort  on  a  closed- 
loop  study. 
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EC-4H  -  Combustion  Instability  Lab 


•  Concentrate  on  near  injector  field 


•  Measure  p’  and  q’  simultaneously  to 
evaluate  the  Rayleigh  criteria  for 
combustion  instability 

•  Start  with  shear  coaxial  jets  -  cold  flow 
heritage 


•  Lab  designed  for  2000  psi  -  about 
double  the  pressure  from  other  labs  in 
the  world 


•  Start  with  current  design  for  acoustic 
drivers  pressure  nodes  and  antinodes 
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Status 
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•  1 500  g  LN2  tank  installed  Sept  28,  201 1 

•  Installed  Class  I  Div  2  outlets  (1 1 5V,  208V,  240V)  -  stripped  floor  for 
02  compatibility,  painted  cell 

•  First  chamber  pieces  arrived  week  of  March  5,  201 2 

•  Fluid  systems  in  full  construction 

•  Data  acquisition  components  starting  to  arrive 


Construction/integration  to  be  finished  by  3Q  2012 
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9  J  is  not  the  only  parameter  to  characterize  mixing  in  shear  coaxial  jets  - 
Area  ratio,  Inner  jet  post  thickness,  and  reduced  densities  matter! 

9  Found  functional  correlations  for  dark  core  lengths 

9  Based  on  four  geometries,  mixing  is  enhanced  with  increasing  J  and  t/DI 
but  less  efficient  with  increasing  Ao/Ai 

9  Through  LES  found  relevant  St  numbers  for  one  geometry  for  different  J 

9  Helical  instabilities  were  observed  for  LAR-thin  and  SAR-thick,  SAR-thin 
geometries  at  large  J  for  baseline  conditions 

9  Geometry  plays  a  role  for  the  response  of  the  jet  to  acoustics,  at  PAN: 

9  LAR-thin:  The  flow  became  less  sensitive  to  applied  PAN  forcing  with 
increasing  J 

9  SAR-thick:  Recirculation  zone  dampens  the  effect  of  PAN  forcing  at  lower 
J,  but  at  higher  J  (>10)  measurable  effect 

9  SAR-thin:  Flow  responds  to  forcing  acoustic  frequency  at  all  J 
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•  Mixing  for  shear  coaxial  jets  has  other  major  variables  other  than  momentum 
flux  ratio  (J) 

—  Ratio  of  inner  jet  temperatures  is  another  important  variable 
—  Geometry  -  Area  ratio  and  lip  thickness  also  affect  mixing 

•  For  LARthick 

—  Found  bending  mode  with  largest  effect  at  velocity  antinodes 

—  The  reduction  on  the  dark  core  length  was  greatest  for  a  medium  J  range 

—  For  near  critical  pressures,  the  collaboration  with  ECP  determined  relevant  St  for  our  injector 
configuration  and  was  able  to  capture  qualitative  behavior  of  natural  and  excited  jets 

•  ForSARthin 

—  Did  not  see  bending  mode  for  conditions  studied 

—  Saw  vortex  roll-up  and  puffing  occurring  over  entire  J  range  (0.09-21)  tested  - 
—  PAN  forcing  produced  symmetric  flow  structures  regardless  of  J 
—  Spectral  plots  showed  strong  response  to  PAN  forcing  at  low  and  high  J 

•  For  LARthin 

—  Sees  both  bending  and  vortex  roll-up  modes  depending  on  the  acoustic  frequency 

—  PAN  forcing  at  low  J  produced  symmetric  flow  structures,  while  at  higher  J,  influence  of  forcing 
subsided 


il  magnitude  plots  showed  decreasing  influence  of  PAN  forcing  with  increasing  J 
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New  Combustion  Chamber 
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AFRL  -  Edwards  was  tasked  to 
develop  an  experiment  and  to  lead  a 
joint  experimental  /  modeling  team 
to  study  coupled  flame  holding 
mechanisms. 

LEVERAGE  EXPERIENCE 
GAINED  WITH  COLD  FLOW 
EXPERIMENTS 
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rWtf 


Construction  and  Organization  of  Data  Set 
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•  First,  form  a  row  vector  consisting  of  all  pixel  intensity  values  of  each 
snapshot  image  (with  resolution  of  n  rows  by  m  columns)  in  order  of 
increasing  columns,  then  increasing  rows 


J  = 


1 

m 


n 

m 


i  =  l  /=1  i  =  2  /  =  2  i  =  n  i  =  n 

i  •  i  O  O  O  •  - g 

7  =  1  7  =  m  7  =  1  7  = m  7  =  1  ./  = m 


•  Then,  combine  all  such  row  vectors  for  sequences  of  image  frames 
resulting  in  a  matrix  A  consisting  of  N  rows  by  (M  =  nxm)  columns  of 

infpnsifv  valnec 
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Matrix  Decomposition 


•  Eigenvalue  decomposition  or  singular  value  decomposition  (SVD)  can  be  used 


•  SVD  Subroutine  readily  available  in  MATLAB® 


•  Prior  to  matrix  decomposition,  the  temporal  mean  of  A  was  subtracted  resulting  in  a  matrix  of  intensity 
fluctuations  A,  i.e.,  i 

A  =  A  --VA  for /=  1...N, 

9  9  NT  9  j=  1...M 


•  Application  of  SVD  on  A  gives  two  orthogonal  matrices  U  (NxN)  and  V  (MxM),  and  a  diagonal  matrix  S 
(NxM)  of  singular  values  in  increasing  order  of  magnitude 


A  =  USV  =  QVt 

•  Thus,  a  time-resolved  set  of  images  intensity  fluctuations  A(jc,t)  can  be  represented  as  a 
linear  combination  of  orthonormal  basis  functions  (j)k  such  that 

M 

a  o,o=2A(o&o) 

k- 1 


where  ak(t )  are  time  dependent  orthonormal  amplitude  coefficients  and  ^(x)  are  proper 
orthogonal  modes  of  A. 


•Equivlance:  columns  of  Q  ~  ak(t),  columns  of  V  ~  ^(x) 
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Acoustic  Analysis  for  Injectors 


In  collaboration 
with  Jeff  Muss  and 
Rory  Davis,  Sierra 
Engineering 


Have  more  accurately 
computed  the  acoustic 
modes  for  the  inner  and 
outer  jets  for  constant  and 
linearly  varying 
temperatures  for 
subcritical  and 
supercritical  pressures 
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Acoustic  Analysis  for  Injectors 


In  collaboration  with  Jeff  Muss  and  Rory  Davis,  Sierra  Engineering 


HJ5TTB-0 

shh  —mas 
shx  -a<mo 


Have  more  accurately 
computed  the  acoustic 
modes  for  the  inner  and 
outer  jets  for  constant  and 
linearly  varying 
temperatures  for 
subcritical  and 
supercritical  pressures 


STEP-l 

SUB  -S 

FRE0-'3159 

['RES  [AVCr> 

H.ETfB-0 

SHH  ”-39087 

SMX  -18&9JJ 
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Injector  IV:  Constant  J ,  Varying  R 


J  =  2.0 

J  =2.0 

rrt 

II 

R  =  3.5 

L!D1  =  17.9 

L!D1  =  12.1 

J  =2.0 

J  =2.1 

• 

^r 

II 

R  =4.3 

L/D1  =  11.6 

L!Dl  =  10.1 
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Collaboration  with  ECP:  Grid  and  Mesh  CENTRALE 

PARIS 
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Boundaries  with  modulated  normal  velocity 


3D  visualization  of  reservoir 


i 

=sin(2rr  f  t ) 


,  Vj=sin(2n  f  t  +  cp) 

Longitudinal  cut  of  domain  with  the  BC's 


Mesh  detail  near  the  injector 


Solver:  AVBP  -  state  of  the  art  LES  code 

*  With  real  fluid  properties  to  tackle  supercritical  fluids 

Mesh:  2  100  000  nodes/10  000  000  tetrahedra 

*  Highly  refined  near  the  injector  (0.032  mm  on  a 
distance  of  10  inner  jet  diameters) 

CPU  hours  on  Europe  Supercomputer 
Center:  100,000 
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Effect  of  Temperature  on  L/D  for  a  given  J 


J  =  9.4  R=9.9 

109K/203K 

Ui=0.93/9.2m/s 

SAR_thinLip 


j  =  9.4  R=9.9 


128K/192K 

Ui=6.6  /36.2m/s 
LARthickLip 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Flow  Structure  Comparison 
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Dark-Core  Length 
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Total  or  curved 
ength  (Lt)  = 

average  of  left 
and  right  sides 


P  =  1.50  MPa  ;  Pr=  0.44;  VR  =  7.50;  J  =  2.64;  acoustic  field  on 


•  Threshold  Images  based  on  Otsu’s  method  (  n.  otsu,  "A Threshold  selection 
Method  from  Gray-Level  Histograms, “  IEEE  Transactions  on  Systems,  Man,  and  Cybernetics,  vol. 
9,  no.  1,  pp.  62-66,  1979.) 

•  Accounts  for  variability  from  image  to  image  (including  dl  the 
parameter  by  which  the  jet  is  normalized) 
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Acoustic  Forcing  Results 


CENTRALE 


\ 


Instantaneous  images  of  the  simulation  of  an  acoustic  case  with  the  injector  at  a  pressure  antinode  for  a  J  =  3.0 
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Synchronized  p’  and  images 
Pr=1.03,  J=5.1,  VR=5.9 
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Synchronized  p’  and  images: 
Pr=1.03,  J=5.4,  VR=6.0 
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Synchronized  p’  and  images: 
Pr=1.03,  J=5.4,  VR=6.0 
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Inner  Injector  Recessed  1  D. 

Pr  =  1,04,  J  =  5.9,  VR  =  3.9 
Forcing  Condition:  3050  Hz,  Maximum  Ap 


March  26,  2010 


U o=33.0  m/s,  Ui=8.50  m/s,  To=184K  ,  Ti=129  K 


MOVIE 
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No  Recess 

Pr  =  1.02, J  ~  5*7,  VR  =  3.7 
Forcing  Condition:  3080  Hz,  Maximum  Ap 


June  21, 2010 


Uo=4.74  m/s,  Ui=1.26  m/s,  To=191K  ,  Ti=130  K 
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Typical  FFT  of  Dark  Core  Length 


Frequency  (Hz)  Time  (s) 


Sampling  Frequency=41kHz,  Driving  frequency=3.05kHz, 
Pch=3.5MPa;  Pr=1.03;  VR=2.68,  MR=2.27, 
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^chamber  =  ^  MPa  (Pt  =  0.44) 


J  =  0.17 


VR  = 

^Ppeak-to-peak  ^Pchamber ' 

Wl  outer/  Wl  inner 


J  =  0.089  with  thin  inner  post  injector 


■,«*«?  '  n  oiumoi-  i  yty-an 


Negligible  effect  vs.  violent  destruction  of  the  jet 

r 
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^chamber  =  ^  MPa  (Pt  =  0.44) 


'  l  :  ‘ 

4, 

EFA 

J  =  2.6  with  thick  inner  post  injector 


J  =  2.0  with  thin  inner  post  injector 


Bending  vs.  vortical  structures  with  fine  atomization 


L 
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Subcritical  Large  J 
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^chamber  =  ^  MPa  (Pt  =  0.44) 


J  =  23  with  thick  inner  post  injector 


J  =  18  with  thin  inner  post  injector 


Two  different  mixing  mechanisms  for  similar  large  J  values 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 
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Nearcritical  Large  J 


^chamber  ^.6  MPa  (Pf  1.1) 


J  =  9.3  with  thick  inner  post  injector 


J  =  9.4  with  thin  inner  post  injector 


LJJkCi  I  CFVtCOYin  ■  Oivuef  I  Mwrt* 


& 


Longer  dark  core  lengths  and  hence  visible  effect  of  acoustics 

for  thin  inner  post  geometry 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Supercritical  Moderate  J 


^chamber  ^’0  MPa  (Pr  1.5) 


J  =  2.4  with  thick  inner  post  injector  J  =  2.6  with  thin  inner  post  injector 


LJJkCi  I  CFVtCOYin  ■  SiVUCf  I 


Large  dark  core  lengths  for  the  thin  post  geometry  prevent  us 
from  observing  acoustic  effects,  if  any 
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Crow  and  Champagne,  1971 

—  Single  jet  preferred  mode,  St=fd/U~0.3 
Ko  etal,  1976-1989 

—  Some  of  earliest  detailed  description  of  near  field  mixing  for  coaxial  jets 
Boldman  et  al,  1975 


Entrommeni  reqon 


Secorriory 

regioi 


—  Experimental  and  theoretical  analysis  for  mixing  of  two  streams  with 
different  velocities  -  points  out  different  vortex  interactions  Kwan  and  Ko>  J 


Sound  and  Vibration,  48  (2),  1! 


Gutmark  and  Ho,  1983 

—  Collects  previous  results  on  jet  preferred  mode,  St  has  a  range  from 
~0.24-0.64 
Dahm  et  al,  1992 


—  Points  out  importance  of  layer  thickness  and  velocity  defect  on  shear  layer 
seminal  pictures  of  different  instabilities 


Villermeaux  1998 

—  Inner  jet  core  length  expression,  LID^=6IJ0  5 
Richecoeur  et  al  (Candel’s  group),  2006 


Dahm  et  al,  JFM  (241)  1992 


—  Forced  transverse  acoustic  excitation  of  flames  | 


Richecoeur  et  al,  JPP  (22)  No  4,  2006 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 
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Subcritical  conditions 

•  Acoustics  undulate 

inner  jet  like  in 
Heidmann’s 
experiments 
consistent  with 
being  driven  by  an 
imposed  velocity 
field 

•  Reductions  of  core 
of  more  than  20% 

for  0.2  <  J  <10 
Near  and  supercritical 
conditions 

•  Less  undulation  on 
inner  jet  but  same 

mode 

•  Not  a  clear  effect  of 

phase  on  the 
reduction  of  the 
dark  core  length 
•  Max  reduction  of 
inner  jet  occurred  at 
1<J<5 


VR  =  2.0  VR  =  5.9 

Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 
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Outer  Jet  Pressure 


Pr  =  1,03,  J  =  5,1,  VR  =  5.9 
Forcing  Condition:  3000  Hz,  Minimum  Ap 


Feb  18,  2010 


Uo=6.  25  m/s,  Ui=1.06  m/s,  To=165K  ,  Ti=118  K 

MOVIE 
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'  attocube 


chamber 

+  Touterje,  (uPstream) 

+  T.  (upstream) 

inner  jet  ^  1  ' 

Outer  Injector  Post 
- Inner  Injector  Post 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Effect  of  Recess:  SAR_ThinLip;  Pr=0.45,  J~0.09 


Flush  inner  post,  J=0.09,  p7p=0.45% 


Recessed  inner  post,  J=0.089,  p7p=0.60% 


Qualitatively  similar  at  very  low  J  values 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Case  I:  Baseline  Flow 
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Mode  1 


Snapshot  Average 


Mode  2  Mode  3  Mode  4 


Mode  5 


Mode  6 
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Mode  1  Mode  2  Mode  3 


Mode  4 


Mode  5  Mode  6 
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LAR_ThinLip 


r 


SAR_ThinLip 


Similar  lip  0! 
thickness/I.D.xSf' 

$5 


LAR_ThickLip 


1 
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f  l  (  rz 

l  l  Vi 
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r . 
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D1 
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t=(D2- 

D1)/2 

Ao/Ai 

0.51 
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1.6 

1.47 

3.96 

4.70 

6.35 

0.84 

2.9 

0.70 

0.89 

2.44 
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Two  recesses 


SAR,  LAR  ->  Small,  Large  Area  Ratio 
ThickLip,  ThinLip  ->  Post  lip  thickness 


LiJkci  I  oiKovm  ■  eirtief  «  muvia  / 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


Chronological  progression  (only  coaxial 
results  are  summarized  in  what  follows) 
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*  Single  jets,  no  coaxial  flow 

—  Davis  et.  al.  (Ph.D.  thesis)  -  single  jets,  no  coaxial  flow 

*  Coaxial  jets 

—  Davis  et  al.  (Ph.D.  thesis)  -  LAR  thick 
—  Leyva  et.al.  -  LAR_thick 

—  Rodriguez  et.  al.  (Ph.D.  thesis)  -  LAR_thick,  SAR_thin 

—  Graham  et.  al.  -  SAR_thin,  two  recesses 

—  Teshome  et.  al.  (Ph.D.  thesis,  expected  April  2012)  - 
complete  all  four  geometries 

9  Also  complete  modal  analysis  of  earlier  geometries 

*  Future:  Combusting  coaxial  jets 

—  Wegener  et.  al.  (Ph.D.  thesis)  (in  process) 

—  Article  219  facility  funds 
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Baseline  Dark-Core  Lengths 
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Image  interpretation  key 
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pressure  =  fixed 
J  (p0u02/ PjUj2)  =  fixed 


Acoustics  on 
0° 

PAN,  VN 

“Pressure  coupled’ 


Subsequent  images: 
acoustics  on; 

C 

Phase  increases  by  45 


Baseline: 

Acoustics 

off 


PN  -  pressure  node  -  Min 
PAN  -  pressure  antinode  -  Max 
VN  -  velocity  node 
VAN  -  velocity  antinode 
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180° 

PN,  VAN 
Largest  difference 

O 

expected  from  0 
“Velocity  coupled” 


360° 

should  look 
similar  to  0° 


Acoustic  Analysis  for  Injectors 


In  collaboration 
with  Jeff  Muss  and 
Rory  Davis,  Sierra 
Engineering 


Have  more  accurately 
computed  the  acoustic 
modes  for  the  inner  and 
outer  jets  for  constant  and 
linearly  varying 
temperatures  for 
subcritical  and 
supercritical  pressures 
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